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Abstract. This paper presents the results of a rheological study of thermotropic nematic colloids aggregated 
into cellular structures. Small sterically stabilised PMMA particles dispersed in a liquid crystal matrix 
densely pack on cell interfaces, but reversibly mix with the matrix when the system is heated above T n i. 
We obtain a remarkably high elastic modulus, G' > 10 s Pa, which is a nearly linear function of particle 
concentration. A characteristic yield stress is required to disrupt the continuity of cellular structure and 
liquify the response. The colloid aggregation in a "poor nematic" MBBA has the same cellular morphology 
as in the "good nematic" 5CB, but the elastic strength is at least an order of magnitude lower. These 
findings are supported by theoretical arguments based on the high surface tension interfaces of a foam-like 
cellular system, taking into account the local melting of nematic liquid and the depletion locking of packed 
particles on interfaces. 

PACS. 61.30.-v Liquid crystals. - 82.70.-y Disperse systems. - 81.40. Jj Elasticity and anelasticity, stress- 
strain relations. 

1 Introduction 

In the preceding paper (jj we have described the struc- 
ture and aggregation mechanisms of thermotropic nematic 
colloid with a new cellular morphology. We argued that, 
in some circumstances, the continuous phase separation 
takes place, driven by the mean field energy of orien- 
tational nematic order in the matrix. This mechanism, 
and the resulting cellular structure with colloid particles 
densely packed on thin interfaces, is different from a num- 
ber of previous observations ^,^,0,0]. In those cases par- 
ticles always aggregated into dense 3-dimensional floes, 
often supported by line or point disclinations. 

We showed that the behaviour of a liquid crystal col- 
loid is determined by the dimensionless parameter WR/K, 
where W is the nematic anchoring energy on the particle 
surface, K the Frank elastic constant]] and R the parti- 
cle radius. A spherical colloid particle in a nematic liquid 
crystal creates a topological defect in the director field, if 
WR/K 3> 1. When this parameter is small (which is our 
case, with PMMA particles of R ~ 150 nm), the parti- 
cle introduces only a small perturbation into the nematic 
matrix. In this case the long range interaction forces be- 
tween particles are sufficiently weak so that the aggrega- 
tion times are long compared with the alternative mech- 
anism of continuous phase separation. Direct observation 
in the bulk of our samples with the confocal microscopy 

For typical thermotropic nematics and homeotropic direc- 
tor anchoring used in our work, these material constants take 
the values W ~ 10~ 6 J/m 2 and K ~ 10" n J/m, e.g. [§. 



technique have confirmed the open-cell morphology and 
the estimate for the cell size of the structure, A ~ (nR)/<l>, 
where nR is the wall thickness and <P the average particle 
concentration, (experimentally we find n ~ 20 — 30, jj)). 
Recently Tanaka et al. jTj] studied phase separation in mix- 
ture of a 5CB nematic and a surfactant, which forms very 
small micelles. They reported some provocative results, 
also quite different from ours. We could speculate that, in 
their case, the colloid parameter WR/K is so small that 
the mechanism of phase separation described in jj) is not 
valid, indeed, nano-micelles might be too small to apply 
any mean-field description of a nematic order. 




Fig. 1. (a) Reflection-mode confocal image of cellular structure 
in the 5CB nematic colloid at 4> = 10% |j| . The white bar shows 
the length scale of 20/xm. (b) The photo of an aggregated 10%- 
5CB colloid, showing its white scattering appearance and the 
evident rigidity. 
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In any case, regardless of which structure the aggre- 
gating liquid crystal colloid adopts, the phase equilibrium 
and kinetics are altered in non-trivial ways due to the un- 
derlying frustrated liquid crystalline order. The cellular 
structure we find in our systems is, clearly, a metastable 
state. One of the main questions we address here is the 
change in rheological behaviour due to locally highly con- 
centrated forces in the aggregated networks, creating high 
barriers to deformation. One may expect a glass-like freez- 
ing of motion and the resulting prolonged stability of the 
phases. 

In contrast to the network supported by topological de- 
fects (disclinations), as in the recent study || of bigger ag- 
gregates and somewhat weak static modulus G' ~ 0.01 Pa, 
the cellular solids in our work possess a storage modulus 
of G' > 10 5 Pa. This solidification was first reported in || 
for a 5CB-PMMA colloid with particles of R » 250 nm. 
In this paper we present further results on the mechan- 
ical properties of nematic colloid, with smaller particles 
and fast quenching rates, showing the storage modulus 
approaching an even higher values, 10 6 Pa, and examining 
the yield and ageing properties of this phenomenon. 

2 Experimental 

The preparation procedures were described in some detail 
in the companion paper A typical thermotropic ne- 
matic liquid crystal 5CB (from Aldrich) has been mixed 
with monodisperse PMMA particles of radius R — 150 nm, 
sterically stabilised by chemically grafted poly-12-hydroxy 
stearic acid chains. Such short radially grafted polymers 
provide a homeotropic director anchoring. In all experi- 
ments we use the relatively small particle volume fraction, 
not more than <P = 15%, determined by weight. Good mix- 
ing of particles was achieved above the nematic-isotropic 
transition and the samples were stored at T ~ 45C (above 
T n ; = 35. 8C in pure 5CB) in a tumbling device to ensure 
that mixtures are homogeneous before any experiment is 
started. 

An alternative colloid system has been prepared with 
another famous thermotropic nematic MBBA (Aldrich). 
We use this system for comparison, as an example of 
a "poor" nematic order in the supporting matrix (the 
MBBA is partially hydrolysed in ambient conditions, with 
effective T ni = 37.3C in pure MBBA), see for detail. 

The mechanical properties have been studied on a Dy- 
namic Stress Rheometer, from Rheometrics Ltd, in the 
parallel plate set-up. A sinusoidal stress of low frequency 
is applied to the sample and the resulting strain is mea- 
sured and interpreted by the software. Unfortunately, this 
device did not have an option to impose an oscillating 
strain of small constant amplitude. This prevented us from 
studying in detail the sample rheology during the nematic- 
isotropic transition and structure formation - the change 
in the modulus is so great (between ~ 0.01 and 10 5 Pa) 
that the constant-stress mode of the rheometer could not 
cope and the settings had to be adjusted below T ni . 

The rheometer was equipped with the purpose-made 
environmental stage allowing the rapid cooling of the sam- 



ple (at a rate 35°/min). The experimental procedure was 
as following: 

1) The sample was heated to ~ 50C in an ultrasonic bath 
(to ensure homogeneous particle mixing) and then trans- 
ferred, via a pre-heated syringe, to the pre-heated rheome- 
ter plate. The upper plate was then lowered slowly, to the 
gap ~ 0.5 mm. All through this manipulation, the tem- 
perature of the sample was remaining above its T n ;. 

2) The environmental chamber was then closed and the 
temperature stabilised at 50C to return the sample to its 
isotropic mixed state. The oscillating movement of plates 
was started at this moment, with a low constant value of 
stress ~ 0.3 Pa since the isotropic suspension is a liquid 
with low viscosity, and a constant frequency of 1 Hz. Sam- 
ples were allowed to equilibrate in this regime for at least 
10 min. 

3) The rheometer run was then stopped and the temper- 
ature quickly dropped to the working level, every time 15 
degrees below T ni (e.g. to 20C for <P = 5% and 18C for 
<P = 15% in 5CB). The rheometer was restarted at ex- 
actly 100 s after the working temperature was reached. A 
higher value of oscillating stress amplitude, ^3 — 30 Pa, 
was required to acquire reliable data on the resulting rigid 
cellular solid. 

4) To study the yield stress behaviour, the same proce- 
dure of sample loading and quenching was followed. The 
rheometer was then set to perform a low frequency stress- 
sweep test, measuring the response of the system on in- 
creasing shear stress - thus reaching the point of breaking 
and yield of the cellular structure. 



3 Results and discussion 

The first observation that one makes on cooling the ho- 
mogeneously mixed colloid suspension below the nematic- 
isotropic transition is the rapid optical change. Above T n j 
the suspension with a low particle concentration is rel- 
atively transparent. Below the clearing point the colloid 
becomes completely opaque, white in the case of 5CB and 
marginally yellow for MBBA. Clearly, a strong multiple 
scattering of light takes place in the material - signifi- 
cantly greater than in a pure nematic liquid crystal, which 
is turbid due to the director thermal fluctuations . Such 
an opaqueness is usually a signature of randomly quenched 
disorder in the nematic director field and is found, for in- 
stance, in nematic silica-gels Q and polydomain nematic 
elastomers [[l0| . Our study of structure of the aggregated 
nematic colloid below T n i in the companion paper p]] con- 
firms this view: the director is anchored on the walls of 
the cellular structure, see fig. |l|(a). 

Secondly, we find the remarkable rigidity of the result- 
ing material. Even before the detailed rheological study, 
one notices that the opaque aggregated colloid is solid. It 
evidently supports its own weight, fig. |l](b), one can cut 
it with a knife, imprint shapes and handle it quite ro- 
bustly. An important unresolved question is of the ageing 
of such a system. In general, the prepared shape of the 
aggregated cellular solid is preserved for a long time, if 
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the temperature is kept constantly deep below T n i. How- 
ever, we observe that some amount of liquid (apparently, 
the pure liquid crystal) leaks from the samples over a pe- 
riod of days and weeks, the effect being much stronger in 
MBBA colloids than in 5CB ones, and at smaller particle 
concentrations.^ In spite of this, the integrity of the sam- 
ple solid shape seemed to be preserved. In this paper we 
do not discuss effects of ageing and leave it for another, 
more detailed research work. 

A third general observation, made both during the un- 
controlled cooling of the colloid mixture on the lab bench 
and in the environmental chamber of the rheometer, is 
that the rate of cooling has a strong effect on the resulting 
rigidity. Fast quenching, at least in the initial period af- 
ter reaching the working temperature, makes the material 
much more rigid and the cellular structure more regular. 
The latter fact is obtained by comparing the confocal mi- 
croscope images of cellular morphology of the same sam- 
ples cooled on the bench and those taken from the rheome- 
ter after a fast quenching and parallel-plate rheometer run. 

Finally, turning to the detailed study of the mechan- 
ical properties of aggregated cellular solids, we find that 
another crucial effect is the actual value of the working 
temperature: the storage modulus G' depends on it very 
strongly. This introduces a slight ambiguity in comparison 
of the data between colloids of different particle concen- 
tration, since they have their T n j different, decreasing with 
the average volume fraction of a colloid mixture <!> M . We 
chose to work at a fixed difference below T ni , but it is 
not clear that the resulting difference in absolute values 
of working T does not introduce a shift in G' . Due to this 
problem, one cannot be certain that the plots of G' against 
colloid concentration, figs. ||(b) and ||(b), are completely 
consistent and one should interpret them with a degree of 
caution. 



3.1 Storage modulus 

After following the same thermal history for all the sam- 
ples, that is a rapid quench to 15° below their respec- 
tive T n j and a 100 s delay after this temperature was first 
reached, the complex modulus G(uj) of the aggregated cel- 
lular solids has been measured as a function of time. We 
are particularly interested in the storage modulus which, 
at low frequency such as uj — 1 Hz, is a representative fea- 
ture of the equilibrium clastic response. Of course, above 
T n i, for a homogeneously mixed low-concentration colloid 
suspension, this storage modulus is essentially zero. 

The resulting values of storage modulus G' are remark- 
ably high for all materials and compositions studied - this 
solidification of a nematic colloid is the main result of this 
paper, illustrated in figs. ^ and |3|. The accurate measure- 
ment of this modulus meets some unavoidable experimen- 
tal difficulties: for a large G' , the stress applied to the 



2 In fact, there is quite a difference between the 5CB and 
MBBA colloids. The MBBA-material is much less robust, 
breaks, leaks liquid and does not have such a long "shelf-life" 
as the 5CB one. 
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Fig. 2. The evolution with time of the storage modulus G' 
at constant frequency ui = 1 Hz for the 5CB colloid at AT ~ 
15° below Tni, for particle concentrations # = 5%, 7% and 
15%, graph (a), curves with increasing saturation. Graph (b) 
shows the dependence of the saturation value G' on the particle 
concentration. A linear fit gives, G' ~ (2 x 10 6 ) <2>Pa. 
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Fig. 3. The evolution with time of the storage modulus G' at 
constant frequency u> = 1 Hz for the MBBA colloid at AT ~ 
15° below T n i, for particle concentrations = 5%, 7%, 10% 
and 15%, graph (a), curves with increasing saturation. Graph 
(b) shows the dependence of the saturation value G' on the 
particle concentration. An attempted linear fit gives G' ~ 
(1.8 x 10 5 ) [<P - 0.045] Pa, 



rheometer plates has to be sufficiently high to produce 
a measurable strain. However, too high a stress results 
in the slipping of the plates, which become lubricated by 
the leaking nematic liquid, or the complete breaking of 
the cellular structure (the yield effect, discussed below). 
The compromise between these two trends leads to rather 
high experimental noise, especially in systems with higher 
particle concentration where the shear strain tends to be 
lower for an allowed range of applied stress. 

The results for G' in 5CB and MBBA aggregated col- 
loids show small variation of the modulus with time, at 
least during the first several minutes or hours. There is 
also a marked difference in values between the 5CB col- 
loid and the "poor" nematic material based on partially 
hydrolysed MBBA. The elastic modulus is more than an 
order of magnitude lower in the MBBA colloid. 

Formation of a rigid cellular solid proven, the final 
level of G' is related to the initial particle concentration. 
Figures |](b) and ||(b) give such a dependence. The data 
points are somewhat scattered and are affected by a high 
noise, or experimental uncertainty. However, the linear de- 
pendence on concentration is a plausible conclusion for 
both 5CB and MBBA colloid aggregates. A striking dif- 
ference between the two systems, apart from the magni- 
tude of G', is the apparent low-concentration threshold at 
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<P ~ 4.5% for the cellular aggregation of the MBBA col- 
loid. This observation in the "poor" nematic model system 
is consistent with the phase diagram in fig. 9 of the com- 
panion paper (l). There, examining the case of smaller 
strength of nematic field, one finds a noticeable region of 
miscibility at low concentrations, below T n j. In contrast, 
for the Landau parameters characteristic of a "strongly 
nematic" 5CB, this region is too narrow to be practically 
observable - accordingly, the G'(<P) dependence for 5CB 
suggests that it aggregates from nearly zero <P. 

At this point one may remark that the rheometric data 
obtained at a slow cooling rate [|| have shown a much 
lower values of storage modulus G' . Although there is a 
difference in colloid particle sizes (the authors of Q used 
R ss 250 nm, which should affect the aggregation kinetics 
and morphology), we also find that the rate of system 
quenching and the time the aggregating colloid spends at 
rest at the working temperature have a strong effect on 
the resulting strength of cellular structure. 



3.2 Stability and Yield stress 

The aggregated cellular structures of nematic colloids were 
very stable over the period of about an hour, as confirmed 
both by visual observation, constancy of confocal images 
and the little variation of measured G' with time. Over 
longer periods, one observes a "leaking" of pure liquid 
crystal from the rather rigid and completely opaque cel- 
lular aggregate. We did not study this effect in detail, but 
one may speculate that the open-cell morphology allows 
the nematic liquid to gradually find its way out of the 
over-constrained cellular environment. 

However, from the moment of formation one observes 
a systematic step-like variation in the modulus, especially 
visible for 15% 5CB system in fig. |^(a). In this particular 
material the G' is so high that we had to apply a relatively 
high constant stress to obtain any measurable deformation 
response. One might argue that the the steps between two 
distinct values of 2.6 and 3.5xl0 5 Pa are due to a partial 
breaking and reforming of the cellular structure.^ 

A different type of rheological experiment, a stress 
sweep at constant low frequency u> = 1 Hz, allows to in- 
vestigate the yield stress behaviour. In figs. [|(a) and ||(a) 
we plot the values of shear strain in response to the grad- 
ually increasing stress between the rheometer plates. The 
slope of the initial linear increase gives the value of G' . 
However, at higher stress the material yields. We identify 
the corresponding value of stress as the yield-stress a y , 
which is a function of colloid concentration. Once again, 
one finds a difference of more than an order of magnitude 
between the 5CB and MBBA materials. 

Plotting the concentration dependence a y figs. [|(b) 
and |(b), we again find that the MBBA colloid shows 
a concentration threshold at <P ~ 4.5% [compare with 

3 An alternative possibility of rheometer plate slipping was 
also frequently observed, but this always resulted in a rapid 
monotonic decline of measured G' to nearly zero - not the 
repeated steps as in fig. |^(a). 
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Fig. 4. The yield stress effect: (a) Strain-stress curves at a 
frequency ui = 1 Hz for the 5CB colloid at AT « 15° below 
T n t, for particle concentrations $ — 3%, 10% and 15%. (b) 
The yield stress vs particle concentration; solid line shows the 
fit Oy w (5400 <£ 1/2 )Pa. 
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Fig. 5. The yield stress effect: (a) Strain-stress curves at a 
frequency w = lHz for the MBBA colloid at AT w 15° below 
Tni, for particle concentrations $ = 5%, 7% and 15%. (b) The 
yield stress vs particle concentration; solid line shows the fit 
a y w (168 [0-O.O45] 1/2 )Pa. 

fig. ||(b)]. A good single-parameter fit to a square- root de- 
pendence cr y oc ^ 1 / 2 is different from the scaling predicted 
by a crude theoretical model in Appendix. However, the 
magnitude of yield stress is in the correct range, appro- 
priately much greater for 5CB colloids where the strong 
nematic field should generate a greater pressure on inter- 
faces and thus preserve the particle dense packing. 



4 Elastic response of cellular superstructure 

First of all, one needs to appreciate the very large elas- 
tic modulus, which usually tells about a significant energy 
density stored in the aggregated nematic colloid. To illus- 
trate the values, let us compare the characteristic energy 
scales involved. Assume, for the purpose of estimate, a 
5% volume fraction of colloid particles with 150 nm ra- 
dius. The director field around each particle contains the 
elastic energy ~ 10 KR in the case of strong anchoring, 
or ~ 0.2 W 2 R?/K in the case of weak anchoring, [ ETj ■ 
Therefore, even the energetically most unfavourable case 
of separate, equidistant particles (separation d ~ 4i? at 
5%) results in the energy density ~ 0.16 K/R 2 ~ 60 Pa. 
In reality, particles would always aggregate in some way, 
significantly reducing this estimate (e.g. particles lumped 
in the nodes of disclination network in |5J are separated by 
d ~ 100 fim, giving the energy density of just ~ 0.01 Pa). 
Also, considering that the particles of this size will have 
the colloid parameter WR/K <C 1, it is far more likely 
that the case of weak anchoring would occur and even the 
separate equidistantly distributed colloid would have the 



V.J. Anderson el ah: Cellular solid liquid crystal colloids: 2. Mechanical properties 



5 



elastic energy density of only ~ 3 x 10~ 3 W 2 / K ~ 10~ 4 Pa. 
None of these values can even crudely account for the large 
values of elastic modulus G' observed in nematic colloids. 

Of course, it is not the energy density itself, but the 
rate of its change - a slope of energy variation with de- 
formation, that determines the modulus. We now try to 
outline the physical origin of such variation. It seems clear 
that the dramatic enhancement of mechanical properties 
is due to the cellular structure of phase-separated nematic 
colloids. Particles, densely packed and solidified on thin in- 
terfaces, produce very high surface tension 7. There are 
several processes contributing to this surface tension. 
(1) Let us first consider the local melting of the nematic 
phase. The gaps between closely packed particles are of 
order or less than £ ~ 10 nm, the characteristic nematic 
correlation length. The nematic liquid in these small gaps 
is melted, cf. fig. ||(a) and the phase diagram in ref. 1 . 
The thermodynamic energy density of a melted nematic 
at T < T n j, AF n is obtained from the mean- field energy 
part of eq. (7) of B; qualitatively 



AF U 



(A a T*) 2 1 - r 
AC t 



10 7 (1 -t) J/m 3 



(1) 



with r = T/T*. The last estimate in ([j]) is based on the 
Landau parameters given in JD. The (1 — T/T*) factor, 
of course, is only valid near the nematic transition, where 
the variation of order parameter Q*(T) is significant. The 
effective surface tension of such an interface is a result of 
increasing the relative volume occupied by the dense phase 
(assuming the constant interface thickness nR) with the 
melted nematic in the gaps between particles: 



7n » (1 - 4>c) (nR)AF n , 



(2) 



where (1 — <f> c ) measures the amount of volume left for 
the nematic liquid in the cell wall {<j) c being the particle 
close-packing volume fraction) . 

(2) The imbalance in energy density between the pure ne- 
matic inside cells on both sides of isotropic interfaces and 
the isotropic melt inside these densely packed interfaces 
acts as an effective pressure of a magnitude ~ AF n . This 
effect would compress the interfaces, trying to reduce the 
volume occupied by the unfavourable melted state. This 
situation is analogous to the process of packing of gran- 
ular matter by gravity or external pressure. Taking our 
monodisperse spherical particles, the results are that the 
maximal particle volume fraction is reached at random 
close packing <f> c « 0.64, when there are on average the 
maximal possible number of contacts, N c w 8.6, between 
the spheres. The contacts create a depletion gap in the free 
energy, see Appendix and fig. ||(b). The resulting cusped 
minimum is very narrow because the dependence of pack- 
ing concentration on the number of contacts is very sharp 
|]l2|,[l3). These, and other theories of granular matter pro- 
vide a number of advanced estimates, some valid over a 
large range of concentrations well below <p c . One may take 
a view that a crude linear estimate N c ~ N ma . x + f3(4> — <p c ) 
can be obtained simply by considering that N c — * al- 
ready at (j) ~ 0-5- This gives the (lower bound) estimate 
for the slope: (3 ~ 100. 




Fig. 6. (a) Model of cell border close-packed with colloid par- 
ticles - the thin interface between clean nematic domains. The 
interface thickness is related to the particle size, d = nR, with 
n ~ 4 in the sketch. The remaining space in the interface is 
filled with the matrix liquid in its isotropic state, (b) The plot 
of free energy F{cj>) (see JlJ] for detail) showing the difference 
in local nematic mean field energy density AF n . The plot also 
shows the cusped minimum of depth AF& resulting from the 
depletion locking between particles in contact (see text and 
Appendix) 



The narrow cusped structure of the free energy at the 
high-</> end of phase diagram provides the additional rea- 
son for the effective surface tension of interfaces, 7d, and 
also gives the expression for the yield stress cr y . The inter- 
face tension is estimated by the following argument: Take 
the value for the energy of a single depletion gap eq. (||) 
and use the estimated linear dependence for the number 
of contacts N c (<j)). We can now describe the local free en- 
ergy density increase in response to the change (decrease) 
in local particle volume fraction (<fi — <p c ) 



F„ 



VR 



l3(-f Ra)(/) c 

VR 



0). (3) 



The particle concentration can be related to the local 
change in surface area 8 A, namely <j}—<j) c K. —<fi c 8A/A c , as- 
suming its constant thickness and the "equilibrium" close- 
packing area A c ^. Then the change in the full free energy 
on stretching the interface can be written as 

(F-F-)^^M^M (4) 



VrAc 



R 2 



7d 



f3j a(j) 2 c X 
R 2 



where vr is the particle volume and A = V/A c the cell 
size. In this way, the eq. (||) identifies the effective surface 
tension 7d resulting from a local energy increase due to the 
rapid loss of particle contacts on stretching the interface. 

We have a cellular solid system with interfaces under 
tension (rather then, say, elastic walls as in juj), with 
open cells, so that the conservation of cell volume is not 
posing a constraint. The clastic modulus of such a system 

4 Take the current density of interface of constant thickness 
nR to be ^ = N/V. Then 



N 



N 



nR(Ac + 8A) nRAc 



<f> c (l - SA/Ac 
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is estimated as G' ~ 7/A, where 7 = j n + 7^ is the total 
effective surface tension of cell walls and the mesh size A ~ 
nR/<P. We thus obtain from the two physically different 
contributions to the surface tension: 



a 



(1 - 4>c) 

[10 6 (1- 



(A T*) 2 (1-t) 

ACt 
-t)<Z> + 10 5 ]Pa 



Pl a<t>l 
R 2 



(5) 



The numerical estimate is obtained for the material pa- 
rameters of 5CB, cf. the Appendix of ref. Q. The ex- 
perimental data of figs. |2|(b) and |^(b) suggest that deep 
below the nematic transition of a classical thermotropic 
nematic material, the dominant factor is the energy of ne- 
matic melting and, thus, the clastic modulus G' is a linear 
function of colloid concentration <f). On the other hand, 
close to T n i, when r — + 1, the modulus should become a 
concentration-independent material constant determined 
by the packing of hard spheres under pressure. 

One of the consequences of the depletion mechanism 
is the effect of the nematic order Q(T), and of small vi- 
brations. This mechanism is based on the particle pack- 
ing approaching the random close packing limit 4> c and 
having a very sharp dependence of average number of de- 
pletion contacts on local concentration. It is known that 
the packing density increases sharply when the sandpile is 
subjected to an external pressure (the role often played by 
the gravity body force) and the vibrations. In our case, the 
nematic mean field creates such an effective pressure. We 
expect, and indeed find in the experimental study of cellu- 
lar structure 0, that a small-amplitude oscillating strain 
in the rheometer makes the cells much more uniform and 
interfaces visibly sharper. 



5 Conclusions 

We have reported results of investigation of mechanical 
properties of new cellular solid aggregates of liquid crys- 
tal colloids quenched below their clearing point. Two main 
conclusions can be made from this study. The magnitude 
of the low-frequency storage modulus, which is compara- 
ble with that of a normal rubber, indicates that the thin 
cell walls possess an unusually high effective surface ten- 
sion. Theoretical arguments suggest that the main mech- 
anism for such a local increase in energy density is the 
local melting of nematic liquid, making G' a linear func- 
tion of average colloid concentration. The second impor- 
tant result of this paper is the yield behaviour of cellu- 
lar macrostructure and the dependence of the yield stress 
on concentration. Qualitative arguments suggest a pos- 
sible mechanism for the observed yield stress. However, 
within the limited accuracy of our experiments, we ob- 
served a concentration dependence er y oc (p 1 / 2 which is 
not accounted for by the crude depletion model. 

An interesting comparison can be made between the 
colloids based on "good" and "poor" nematic matrix. Not 
only the values of elastic modulus and the yield stress are 
much lower in the MBBA system, but we also observe 
a characteristic threshold concentration <P ~ 4.5% below 



which the MBBA colloid does not phase separate into the 
cellular structure. 

We must emphasize a number of experimental difficul- 
ties in the study of mechanical properties of our cellular 
solids. The rheometric results, even for the same batch of 
colloid, are often not reproducible. The reasons are the 
important role of the rate of quenching and, specifically 
for the parallel-plate rheometric experiment, slipping of 
plates lubricated by the nematic liquid. It is possible that 
an alternative approach, for instance examining the me- 
chanical response to an indentation probe rather than the 
sheared plate, will give a better defined experiment. 

Notwithstanding these uncertainties, we have observed 
the unusual phase behaviour and remarkable rheological 
response in the "classical" liquid crystal colloid - a straight- 
forward low-concentration mixture of thermotropic nematic 
with sterically stabilised PMMA particles. The metastable 
cellular solid aggregates are preserved by thermodynamic 
forces provided by the nematic mean field. No doubt, fur- 
ther investigations into the mechanisms of cellular struc- 
ture formation and its stability would make possible a 
controlled, reliable and robust preparation of a wide vari- 
ety of new materials. 

We appreciate valuable discussions with M.E. Cates, S.M. 
Clarke and P.D. Olmsted. This research has been sup- 
ported by EPSRC UK. 



Appendix: 

Depletion gap in densely packed walls 

Another effect that switches on when the particles come 
into close contact on localised interfaces is the depletion 
reduction of their surface energy. When isolated in the 
nematic solvent matrix, each particle's surface has the 
energy of the order 7 D (47ri? 2 ), with 7 Q the main surface 
tension of, e.g. hydrocarbon PHSE and biphenyl oil inter- 
face; qualitatively the value of this tension is of the order 
7o ~ 0.05 J/m 2 [^5) and far exceeds all anisotropic an- 
choring corrections. In the Deryagin approximation, each 
contact between two particles creates a depletion disk and 
reduces this energy by ~ 2~f Ra, where a ~ 1 nm is the 
size of the molecule of suspending liquid. In the randomly 
close-packed state each particle has iV c ~ 8.6 contacts (the 
never-reached absolute maximum would be N c — 12 for 
fee packing) which, for a cell wall of thickness nR, gives 
the depletion gain in the local potential energy density 



AF d ~N c (-Y Ra)j^ 



5 x 10 3 J/ni 



(6) 



This gap would continuously deepen on increasing the 
number of contacts. However, the maximal number N c 
is reached at close packing and cuts this linear decrease, 
thus creating a cusped minimum of the free energy density, 
shown in fig. [|(b). The depth of the minimum (^) gives 
an estimate for the yield stress: when the interface is dis- 
torted such that the particles lose most of their contacts, 
the high resistance to deformation is lost and the cellular 
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structure would no longer be stable. Taken in proportion 
to the total amount of interface in the cellular structure, 
the macroscopic yield stress should be a function of initial 
concentration of the colloid: 

cr y - <P AF d ~ 200 Pa. 

One should appreciate that the estimates here are ex- 
tremely crude, and yet they give not an unreasonable value 

CTy < G 1 . 

References 

1. V.J. Anderson, E.M. Terentjev, S.P. Meeker, J. Crain and 
W.C.K. Poon, the companion paper on "morphology". 

2. P.E. Cladis, M. Kleman and P. Pieranski, C. R. Acad. Sci. 
Paris 273 (1971) 275. 

3. P. Poulin, V.A. Raghunathan, P. Richetti and D. Roux, J. 
Physique II 4 (1994) 1557. 

4. P. Poulin and D.A. Weitz, Phys. Rev. E 57 (1998) 626. 

5. M. Zapotocky, L. Ramos, P. Poulin, T.C. Lubensky and D. 
Weitz, Science 283 (1999) 209. 

6. P.G. de Gennes and J. Prost, Physics of Liquid Crystals 
(Clarendon, Oxford, 1993). 

7. H. Tanaka and J. Yamamoto - to be published (2000) 

8. S.P. Meeker, W.C.K. Poon, J. Crain and E. Terentjev, Phys. 
Rev. E (2000) accepted. 

9. T. Bellini, N.A. Clark, CD. Muzny, et al., Phys. Rev. Lett. 
69 (1992) 788. 

10. E.M. Terentjev, J. Phys. C: Cond. Matter 11 (1999) R239. 

11. O.V. Kuksenok, R.W.Ruhwandl, S.V. Shiyanovskii and 
E.M. Terentjev, Phys. Rev. E 54 (1996) 5198. 

12. S. Torquato, Phys. Rev. Lett. 74 (1995) 2156. 

13. J.X. Liu and T.J. Davies, Powder Metallurgy 40 (1997) 
48. 

14. M.F. Ashby and L.J. Gibson, The mechanical properties of 
cellular solids (Cambridge University Press, 1983). 

15. Treatise on Solid State Chemistry, ed. N.B. Hannay, Vol. 
6B - Surfaces II (Plenum Press, N.Y. - London, 1976). 



